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Abstract 



Ab initio calculations of the electronic energy loss of ions moving in alu- 
minum crystal are presented, within linear-response theory, from a realistic 
description of the one-electron band-structure and a full treatment of the dy- 
namical electronic response of valence electrons. For the evaluation of the 
density-response function we use the random-phase approximation and, also, 
a time-dependent extension of local-density functional theory. We evaluate 
both position-dependent and random stopping powers, for a wide range of 
projectile velocities. Our results indicate that at low velocities band struc- 
ture effects slightly enhance the stopping power. At velocities just above the 
threshold velocity for plasmon excitation, the stopping power of the real solid 
is found to be smaller than that of jellium electrons, corrections being of about 
10%. This reduction can be understood from sum rule arguments. 
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I. INTRODUCTION 



The stopping power for charged particles penetrating a sohd has been the topic of con- 
siderable theoretical and experimental interest, since the beginning of this centuryEHi. The 
electronic stopping power due to collisions with valence electrons has been evaluated for 
many years on the basis of a jellium model&i of the target, the electronic states being de- 
scribed by plane waves. However, in a more realistic approach valence electrons move in a 
periodic potential, electronic states are described by Bloch wave functions, and the spec- 
trum of one-electron excitations splits into the so-called energy bands. The impact of band 
structure effects on both plasmon dispersion curve si and dynamical structure factorsHlSl has 
been investigated recently, demonstrating the importance of these effects even in the case of 
free-electron-like metals such as aluminum. 

Among the most recent attempts to fully introduce the electronic band structure in the 
evaluation of the stopping power for low projectile velocities there is, for alkaline metals, 
a one-band calculation0, as well as a calculation based on a linear combination of atomic 
orbitalJii. For arbitrary incident velocities the stopping power can be calculated, within 
linear-response theory, from the knowledge of the dynamical density-response function of 
the targetEHli. Approximate semiempirical treatments of this quantity have been made, 
and stopping powers of siliconlli^ and goldS for channeled ions have been predicted. More 
recently, the low velocity limit has been investigated, on the basis of a static treatment of 
the density-responseiil. First-principles treatments of the full dynamical electronic response 
of various solids have also been performed, and preliminary ab initio evaluations of the 
stopping power of real solids have been presentedil^. 

In this paper we investigate, within linear-response theory, the valence electronic energy 
loss of ions moving through aluminum crystal. Aluminum is well known to be a typical 
jellium-like metal with a well-defined excitation spectrum, and it represents, therefore, an 
appropriate benchmark for ab initio evaluations of the electronic stopping power of real 
solids. First of all, in section II we describe our full treatment of the wave-vector and fre- 
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quency dependent electronic response of valence electrons, based on a realistic description of 
the one-electron band structure and first-principles pseudopotential theory. For the evalua- 
tion of the density-response function we use the random-phase approximation (RPA)@ and, 
also, a time-dependent extension of local-density functional theory (TDLDA)0. In section 
III, we derive explicit expressions for both random and position-dependent stopping powers, 
from the knowledge of the imaginary part of the projectile self-energy, which is evaluated 
in the so-called GW approximationil. In section IV, numerical calculations of the stopping 
power of valence electrons in aluminum crystal are presented, for a wide range of projectile 
velocities, and we compare our results with the stopping power of a homogeneous electron 
gas with a density equal to that of aluminum. We evaluate, separately, contributions to the 
energy loss coming from the excitation of single electron-hole pairs and plasmons, and we 
interpret our results on the basis of the f-sum-rule for the dynamical structure factor. We 
only consider the electronic response of valence electrons, and the contribution to the elec- 
tronic stopping power coming from the excitation of core electrons is, therefore, not taken 
into account. In section V, our conclusions are presented. 

II. DENSITY-RESPONSE FUNCTION 

The linear density-response function x(r,r',a;) of an electron system is defined by the 
equation 



where p*"'^(r, cu) is the electron density induced by an external potential V'^^'^{r,uj). 

In a self- consistent- field theory, the induced electron density is derived from the response 
function r', u) for non- interacting electrons moving in an effective potential Veff{r,u), 
as follows 




(2.1) 



dh'x^T, r', u) [r"^*(r', u) + V^"<^(r', u)] 



(2.2) 
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where l/*"''^(r, cj) represents the hnear change in Veff{r,u) brought about by the induced 
electron density itself. Since we consider a time-dependent external field, we are in the more 
general scenario of time-dependent density functional theory (DFT)ii, whose theoremsc3 
generalize those of the usual DFT0. 

Within time-dependent DFT, V'^"''^{r, u) consists of the sum of two terms (we use atomic 
units throughout, i. e., rrie = e = h = 1): The Hartree contribution, 

\/A"'(r,^) = / dV^-p*"'^(r',^), (2.3) 
J |r — r I 

which accounts for the average (long-range) effects of the Coulomb interaction between the 
target electrons, and the exchange-correlation (XC) contribution, 

V:f{r,co) = I d^r'^^p-'^(r',.;), (2.4) 

which accounts for the effects of all many-body short-range correlations not included in the 
Hartree approximation. Here, Vxc represents the functional derivative of the XC energy 
functional E. 



-'XC- 



VM = (2.5) 



Introduction of Eqs. (U) and (U) into Eq. and comparison with Eq. (|2.1j ) result 



in the following integral equation for the self-consistent-field density-response function: 

X-(r. r'; = x>. r'; -) + / dr, / dr.xAr. .,,.)V(r. r,, .)x(r.. (2.6) 



where 



V^(r,r» = ^- + ir,,[p] (2.7) 



r — r 



and 

'^-W = , ,^'u°'fl V P.8) 

Sp[r,u;)6p{r',uj) 

In the time-dependent Hartree or random-phase approximation, Ve//(r) consists only of 
the average electrostatic interaction between the electrons, the induced potential is then 



given by the Hartree term exclusively, and K^dr, r') in Eq. (|2.7| ) is taken to be zero. In the 
TDLDA, the kernel K^d^^, r') entering Eq. ( |2.7| ) is replaced by 



ir.,(r,r')^^^ = 5(r-r') 



(2.9) 

P=Po(r) 



dp 

where Krc(p) is the derivative of the XC energy of a homogeneous electron gas of density 
p, which we compute with use of the Perdew and Zunger parametrization0. This so-called 
TDLDA represents an adiabatic extension (a zero frequency K^c is used) to finite frequencies 
of the local density approximation (LDA) for XC. 

For periodic crystals we introduce a Fourier expansion of the density-response function, 

1 BZ 

X{r,r',uj) = E e^^^^''^-^e-(^+^')-r'xG,G'(q,c^), (2.10) 

q G,G' 

where Q represents the normalization volume, the first sum runs over q vectors within the 
first Brillouin zone (BZ), and G and G' are reciprocal lattice vectors. Then, introduction 
of Eq. (|2.10| ) into Eq. ( p. 61) leads to the following matrix equation: 

XG,G'(q,^^) = XG,G'(q)^) + EE^G,G"(q,^)^G",G"'(q,^^)xG"',G'(q,^), (2-ii) 

G" G'" 

where VG",G"'(q, i-^) represent Fourier coefficients of the interaction potential of Eq. (^.71). 
The Fourier coefficients of the density-response function of non-interacting electrons have 
the well-known formil 

/ \ 1 /k,n /k+q,n' 

XG,G'(q>^) = n EE 



X (0k,n I e~^('l+^)-'- 1 0k+q,n' ) (0k+q,n' | c'^^^'''^'' \ 0k,n) , (2.12) 

where the second sum runs over the band structure for each wave vector k in the first BZ, 
/k,n are Fermi factors, and (/)k,n(r;^) are single-particle Bloch states. 



In the RPA, the one-electron Bloch states entering Eq. (|2.12|) are the self-consistent 



eigenfunctions of the one-electron Hartree Hamiltonian, and the coefficients VG",G"'(q, t^) 



entering Eq. (|2.11| ) are the Fourier coefficients of the bare Coulomb interaction. In the 



TDLDA, the one-electron Bloch states entering Eq. ( |2.12|) are the self-consistent LDA 
solutions of the Kohn-Sham equation of DFT, and the coefficients Vg",g"'(c1; ^) entering 
Eq. ( |2.11| ) are the Fourier coefficients of the interaction potential of Eq. ( |2.7| ) with the XC 
kernel of Eq. (^). 

For the evaluation of the one-electron Bloch states, we first expand them in a plane wave 
basis, 



-^X:«M(G)e^(''+«)-, (2.13) 

with a kinetic energy cutoff of 12 Rydbergs, which corresponds to keeping approximately 100 
G- vectors in Eq. ( |2.13|) . Then we solve for the coefficients Mk,n self-consistently, within a full 
description of the electron-ion interaction based on the use of a non-local, norm- conserving 
ionic pseudopotentialii. We subsequently evaluate, from Eq. ( |2.12|) , the Fourier coefficients 
Xg g'(^' which is the most demanding part of the response calculation. Finally, we solve 



the matrix equation ( |2.11| ) for the Fourier coefficients of the interacting response function 
Xg,g'(^! C(j), which we obtain in both RPA and TDLDA. 



III. ELECTRONIC STOPPING POWER 

We consider a particle of charge Zi moving in an inhomogeneous medium. The damping 
rate of the projectile in the state 0o(r) with energy Eq is obtained from the knowledge of 
the imaginary part of the self-energy S(r, r'; Eq), according to: 

7 = -2 1 dr 1 dr'0*(r)ImE(r, r'; i?o)0o(r'). (3.1) 

In the so-called GW approximation, the self-energy is given byS 

i r'^ 

E(r,r';w) = — / du;'G{r,r';u - uj')W{r,r';uj'), (3.2) 

271 J-oo 

where G{r,r';uj) represents the Green's function for the projectile, and W{T,r',uj) is the 
screened interaction: 
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W{r, r'; u) = v{r - r') + J dri J dr2t;(r - ri)x(ri, V2, uj)v{r2 - r'), (3.3) 

v{r — r') being the bare Coulomb interaction and x(r,r',Ci;) the density-response function 
of the medium. For periodic crystals, we introduce the Fourier representation given by Eq. 
( ^.l(j| ) into Eq. ( |3.3[ ), and we find that 

y2 BZ 

Wir,r',uj) = -^J2 E e^^''^''^-'-e-(''+^')-^V(q)[5G,G' + XG,G'(q,^)^G'(q)], (3-4) 



q G,G' 
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where VG{q) = 47r/|q+ G| are the Fourier coefficients of the bare Coulomb interaction. 

Replacing the Green's function entering Eq. ( p.2|) by its zeroth order approximation, we 
find the following expression for the damping rate: 
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7 = / duP^, (3.5) 
Jo 

where 

Pu = -2Y.jdv j dr'0}(r')0S(r)W(r,r',Eo -i5;/)0/(r)0o(r')5[cc'- (i?o -i^;/)], (3.6) 

the sum running over a complete set of final states 0/(r) of energy Ej. 

The quantity P^^ entering Eq. ( |3.5| ) represents the probability for the projectile to transfer 
energy uj to the medium. Consequently, the stopping power of the medium, i. e., the energy- 
loss per unit path length of the projectile, can be obtained as follows: 

dE 1 r°° 

duuP^, (3.7) 



dx V Jo 

where P^^ is given by Eq. 



A. Random stopping power 

In the case of heavy projectiles moving at a constant velocity v with no definite impact 
parameter, the initial and final states can be described by plane waves: 

Mr) = -^e'^-^ (3.8) 



and 



0/(r) 



^i(qo-q)-r 



(3.9) 



where qo and q represent the initial momentum of the projectile and the momentum transfer, 
respectively. 

For heavy projectiles recoil can be neglected, i.e., 

Eo-Ef = q-v, (3.10) 

and after introduction of Eqs. ( |3.4| ), ( |3.8| ), ( p.9| ) and ( |3.10| ) into Eq. ( |3.6| ) we find the 
following expression for the so-called random stopping power: 



dE_ 

dx 



- random 



where eQ q, (q, uj) is the inverse dielectric matrix in momentum space: 

^G,G'i^^^) = ^G,G' + ] — ^XG,G'(q,^), 

|q + 

the density-response matrix XG,G'(q>^) being given by Eq. ( p.ll| ). 

The symmetry of the one-particle Bloch states results in the following identity: 

(^G^G'iSq, uj) = e5-iG,5-iG'(q,t^), 



(3.11) 



(3.12) 



(3.13) 



S representing a point group symmetry operation in the periodic crystal. As a consequence. 



the stopping power of Eq. ( 3.11 ) can be evaluated from the knowledge of the dielectric 
matrix for wave vectors q lying in the irreducible element of the Brillouin zone (IBZ): 

_ o„ IBZ 



dx 



= E E E tI^^^i-^5^^g,5-g [q, (^q + O ■ v] , (3.i4) 

random q 5 G Pq+*^l 



the second sum in this equation running over the symmetry operations generating the wave 
vectors in the star of each q. 

If the diagonal elements of the inverse dielectric matrix entering Eq. ( |3.11| ) are replaced 
by the inverse dielectric function of a homogeneous electron gaJl, i. e.. 



^homog 



(Iq+GU), 



(3.15) 



Eq. ( ^.111 ) exactly coincides with the well-known formula for the stopping power of a homo- 
geneous electron ga^. 
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B. Position-dependent stopping power 



In the case of heavy projectiles moving with constant velocity v on a definite trajectory 
at a given impact parameter b, the initial and final states can be described in terms of 
plane waves in the direction of motion and a 6 function in the transverse direction!!. Then, 
introduction of these states, 0o(r) and 0/(r), and the screened interaction of Eq. ( |3.4| ) into 
Eq. ( P^ ) gives, after neglecting the projectile recoil, the following result for the position- 
dependent stopping power: 



dE 
dx 



(3.16) 



the sum ^ ' being restricted to those reciprocal lattice vectors that are perpendicular to 

K 

the velocity of the projectile, i. e., K • v = 0. eQ^Q,(q, cj) is the inverse dielectric matrix of 
Eq. ( p.l2| ) and S represents, as in Eq. ( p.l4| ), a symmetry operation of the point group of 



the periodic crystal. 

The most important contribution to the position-dependent stopping power of Eq. ( p.l6| ) 



is provided by the term K = 0, the magnitude of the other terms depending on the direction 
of the velocity. For those directions for which the condition of K ■ v = is never satisfied we 
have the random stopping power of Eq. (|3.14|) , and for a few highly symmetric or channeling 
directions non-negligible corrections to the random result are found, thus exhibiting the 
anisotropy of the position-dependent stopping power. We also note that (a) the average 
over impact parameters of the position-dependent stopping power of Eq. ( p.l6|) along any 



given channel has the same value as the random stopping power of Eq. ( p.l4| ), and (b) as 
long as the diagonal elements of the so-called dynamical structure factor, — 2 ImxG,G(ct; i^), 
were isotropic, there would be no dependence of the random stopping power on the direction 
of V. 



9 



IV. RESULTS AND DISCUSSION 



The input of our calculation of both random and position-dependent stopping powers 
of Eqs. (3.14) and (3.16) is the interacting response matrix Xg,g'(^)^) of Eq. (2.11), 
which we have solved for Al crystal. The number of bands required in the evaluation of the 
polarizability XQQ,(q, cj) of Eq. (2.12) depends on the value of the frequency oj. As the 
maximum energy uj transferred by the moving projectile to the target is proportional to the 
velocity, only a few bands are required in the case of slow projectiles. Indeed, it has been 
foundii that at velocities below v = O.la.u. just two bands are required to account for the 
energy-loss. At larger velocities excitations to higher bands become possible, and we find 
that convergence is achieved, for v < 2a.u., if one considers on the order of 30 bands. The 
results presented below have been found to be well-converged for all velocities under study, 
and they all have been performed with the use of 60 bands. 

The sampling of the BZ required for the evaluation of both the polarizability of Eq. 
(2.12) and the stopping powers of Eqs. (3.14) and (3.16) has been performed on 10 x 10 x 10 
Monkhorst-Pack meshesii (47 points in the IBZ) with an energy smearing of the Fermi level 
of 0.25eV, and this sampling has been done in conjunction with a choice of a finite value of 
the damping factor entering Eq. (2.12) of rj = 1.5eV. At velocities of the projectile below 
the velocity for which plasmon excitation becomes possible, convergence is already achieved 
with the use of 4 x 4 x 4 Monkhorst-Pack meshes (8 points in the IBZ). 

A. Random stopping power 

The main ingredient in the calculation of the random stopping power of Eq. (3.14) is the 
so-called dynamical structure factor SG{q^,uj). It is proportional to the diagonal elements of 
the imaginary part of the density-response matrix: 

S{q+G,uj) = -2ImxG,G(q, ^) 
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'^ + ^'-Wc.(q,^), (4.1) 



27r 

where q is taken to be in the first BZ. Thus, couphngs of the wave vector q + G to wave 
vectors q + G' with G 7^ G', which appear as a consequence of the existence of electron 
density variations in real solids, only contribute to the random stopping power through the 
dependence of the diagonal elements of the interacting response matrix Xg,g(cI;|^) on the 
off-diagonal elements of the polarizability Xg g'I^' have found that in the case of Al 

crystal, which does not present strong electron density gradients nor special electron density 
directions (bondings), contributions from these so-called crystalline local field effects are 
within 0.5% of the total random stopping power. Thus, as for the evaluation of the random 
stopping power of Al crystal, the off-diagonal elements of Xg g' (^j ^) '^^^ t)e neglected. 

In the general case of an inhomogeneous system one finds, for the imaginary part of the 
inverse dielectric function, the following form of the f-sum-rule: 

roc 

/ dco'cuImeQ^'^Q, (q, cj) = — 27r^r?,G-G', (4-2) 
Jo 

where ug represent the Fourier components of the density. Note that ^0=0, which equals 
the average density of the system, does not depend on the details of the band structure. 
Thus, it may be argued that band structure effects, which have an impact on both plasmon 
dispersionS and the dynamical structure factorS, may give no corrections to integrated 
quantities as the random stopping power. However, because in the sum over the wave 
vector q in Eq. (3.14) the frequency u in the argument of the structure factor takes values 
from to only |q + G| f , the stopping power of band electrons and that of electrons in an 
effective jellium with the same average electron density will be, in general, different. 

In Fig. 1 we show, as a function of the velocity of the projectile, our full calculation of the 
random stopping power (see Eq. (3.14)) of valence electrons in Al crystal for protons {Zi = 
1), together with the corresponding result for a homogeneous electron gas with an electron 
density parameter equal to that of aluminum (r,, = 2.07)0. In the calculation of the stopping 
power of both band and jellium electrons, the third sum in Eq. (3.14) is extended over 15G 
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vectors of the reciprocal lattice, the magnitude of the maximum momentum transfer q + G 
being 2.9qF {qp is the Fermi momentum Solid and open circles represent our calculated 
random stopping powers of band electrons, as calculated within the RPA and the TDLDA, 
respectively. Solid and dashed lines represent the corresponding stopping powers of jellium 
electrons. All these results have been found, for Al crystal, to be insensitive to the choice 
of the projectile velocity direction. 

In the case of projectiles moving at low velocities, below the threshold velocity for plas- 
mon excitation, the sum over the frequency = (q + G) ■ v in Eq. (3.14) can never be 
replaced by an integration from to oo as in Eq. (|4.2|) , and the stopping power will, in prin- 
ciple, depend on the band structure of the crystal. The existence of interband transitions, 
not present within a jellium model of the crystal, result in a dynamical structure factor 
which is, at low frequencies (where it increases linearly with frequency), slightly enhanced 
with respect to the corresponding jellium result0. As a result, the stopping power of the 
real target is, for projectile velocities smaller than the Fermi velocity and within both RPA 
and TDLDA, a linear function of the velocity and about 7% higher than the stopping power 
of jellium electrons. 

The threshold velocity for which plasmon excitation becomes possible is Vt ~ 1.3a.u. 
in the case of electrons in jellium with = 2.07. In the case of Al crystal this threshold 
velocity is slightly smaller {v ~ 1.2a.u.), as can be easily concluded from an inspection 
of the plasmon dispersion curve for Al crystal of Ref. ^ For projectile velocities larger 
than Vt there are two mechanisms of valence electron excitation: (a) excitation of single 
electron-hole pairs and (b) excitation of plasmons. Figure 2 exhibits, as a function of the 
projectile velocity, separate contributions to the stopping power of both band (solid circles) 
and jellium (solid lines) electrons coming from electron-hole pair excitation (Fig. 2a) and 
plasmon excitation (Fig. 2b) . Although for wave vectors that are smaller than (the critical 
wave vector where the plasmon dispersion enters the electron-hole pair excitation spectrum) 
both mechanisms of valence electron excitation contribute to the energy-loss, contributions 
from losses to electron-hole pair excitations are negligible for q < qc, and Eq. (4.2) leads 
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us to the conclusion that contributions from losses to plasmon excitation are independent 
of the detailed band structure of the crystal. That this is the case is obvious from Fig. 
2b. As for the contribution to the energy-loss coming from the excitation of electron-hole 
pairs, Figs. 1 and 2a show that band structure corrections lower the stopping power of 
electrons in jellium by about 10% at and just above the plasmon threshold velocity; this is a 
consequence of the dynamical structure factor for the real crystal being, within the electron- 
hole pairs continuum at q > Qc and u < {q + G)v, smaller than in the case of a homogeneous 
electron gas. Nevertheless, contributions to the stopping power coming from the excitation 
of plasmons are still smaller than contributions from losses to electron-hole pairs, and there 
is, at high velocities, exact equipartition of the energy-loss, as in the case of a homogeneous 
electron ga Also, at high velocities the sum over the frequency = (q + G) ■ v in Eq. 
(3.14) can be replaced by an integration from to oo and the stopping power of both jellium 
and band electrons is expected to be the same, as illustrated in Figs. 1 and 2a. 

We note that the TDLDA results plotted in Fig. 1 for the stopping power of jellium 
electrons, which describe short-range correlation effects within the LDA, approximately 
reproduce more detailed calculations with either static or dynamic local-field corrections 
included^. The stopping power of both jellium and band electrons is enhanced by about 
20%, at low velocities (exclusive electron-hole pairs domain), with respect to the RPA result, 
as a consequence of short-range correlation effects provoking a reduction in the screening 
within the target. The sum rule of Eq. ( |4.2|) equally applies for both RPA and TDLDA 
dynamical structure factors, and this prevents, therefore, plasmon contributions to the stop- 
ping power from being sensitive to the details of the wave- vector and frequency dependence 
of the response. At large velocities, well above the stopping maximum, all calculations pre- 
sented in Fig. 1 converge to the well-known Bethe formulaS, in which the mean excitation 
potential is replaced by the plasma energy@. 



13 



B. Position-dependent stopping power 



We have carried out, from Eq. (3.16), calculations of the electronic energy-loss versus 
impact parameter, and we have obtained, within the RPA, the results presented in Fig. 3. In 
this figure we represent by solid curves the position-dependent stopping power of Al crystal 
for protons {Zi = 1) moving with v = 0.2a.u. along the (100) (Fig. 3a) and (111) (Fig. 3b) 
directions. In Fig. 3a the impact parameter has been taken to be b = (0, b, 0)ac (oc is the 
lattice constant), with b from to 1, and in Fig. 3b it has been taken to be b = (—1 1 0)ac. 
Contour density plots of the square lattice containing both the projectile trajectory and the 
b vector are displayed in Figs. 4a and 4b for the (100) and (111) directions, respectively. 
The averaged electronic densities along the projectile paths depicted in Fig. 4 are plotted 
in Fig. 5, as a function of the impact parameter, showing for this fee crystal the expected 
periodicity with b = 1/2, also present in the position-dependent stopping power of Fig 3. 

First of all, we note that the existence of small electron density variations in real alu- 
minum results, through the off-diagonal elements of the interacting density-response function 
Xg,g'{^.,^), in non-negligible differences between position-dependent and random stopping 
powers. We have obtained differences up to 20% for projectiles incident in the (100) direction 
(see Fig. 3a), and up to 10% for projectiles moving in the (111) direction (see Fig. 3b). The 
maxima in the stopping power for trajectories along the interstitial regions (6 = 1/4) and 
the minima near the cores (b=0) are associated with corresponding maxima and minima in 
the integrated electronic densities of Fig. 5. 

A comparison with a local density approximation (LDA)ii is also shown in Fig. 3 (dashed 
curves). In this approach, the position-dependent stopping power is obtained as the stopping 
power of a homogeneous electron gas with an electron density equal to the average electron 
density along the projectile path. Although the impact parameter dependence, as obtained 
in this approach, is qualitatively the same as found within our full band structure calcu- 
lation, position-dependent corrections to the random stopping power appear to be largely 
underestimated within the LDA. This is obviously a consequence of the stopping power 
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depending not only on the local density, but also on the electron density gradients further 
away from the projectile trajectory. We have also calculated position-dependent stopping 
powers of Al crystal for various non-relativistic velocities of the incident projectile, and we 
have found that corrections to the random stopping power, as obtained within the LDA, 
are too small for all velocities. These corrections are found, within our full treatment of the 
band structure, to slightly decrease with the velocity of the projectile. In the case projectiles 
moving near the cores {b — 0) along the (100) direction, we have obtained differences be- 
tween position-dependent and random stopping powers of approximately 20% — 12% in the 
velocity range from v = 0.2a.u. to v = la.u.. For projectiles moving near the cores {b = 0) 
along the (111) direction, these differences have been found to be, in the same velocity range, 
of approximately 10% — 8%. 



V. CONCLUSIONS 

We have presented full band structure calculations of both random and position- 
dependent stopping powers of valence electrons in Al crystal, as obtained within linear- 
response theory and with the use of the RPA and the TDLDA. The random stopping power 
has been evaluated for a wide range of projectile velocities. Our results indicate that at 
low velocities of the projectile band structure effects result in the stopping power of real 
aluminum being about 7% higher than that of a homogeneous electron gas with a density 
equal to the average electron density of the real crystal {vg — 2.07). At velocities just above 
the threshold velocity for plasmon excitation, the stopping power of valence electrons in Al 
crystal is found to be smaller than that of jellium electrons, corrections being of about 10%. 
In the high-velocity limit, the stopping power of both jellium and band electrons is found 
to be the same and to coincide with the well-known Bethe formula. 

The various contributions to the random stopping power of Al crystal have been calcu- 
lated separately, showing that band structure effects on the contribution to the stopping 
power coming from plasmon excitation are negligible. This result has been shown to be a 
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consequence of the dynamical structure factor of both jeUium and real aluminum fulfilling 
the same sum rule, as long as the average electron density is the same. Also, there is, in the 
high-velocity limit, exact equipartition for the energy-loss, as in the case of a homogeneous 
electron gas. 

The position-dependent stopping power has been evaluated for projectiles with veloci- 
ties up to the Fermi velocity incident in various high-symmetry directions. We have found 
differences between position-dependent and random stopping powers up to 10% for pro- 
jectiles incident in the (100) direction and up to 20% for projectiles moving in the (111) 
direction. The magnitude of these position-dependent relative corrections to the random 
stopping power has been found to be, within the LDA, largely underestimated, specially at 
high velocities of the projectile. 
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The Bethe formula is, in the high velocity limit, nearly insensitive to differences among the 
plasma frequencies obtained in the various approximations to the response of the target. 

The use of this so-called local density approximation (LDA) should not be confused with 
the use, within DFT, of the LDA for the evaluation of T4c and the kernel K^c of Eqs. (2.5) 
and (2.8). 
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FIGURES 

FIG. 1. Full band structure calculation, from Eq. (3.14), of the random stopping power of 
valence electrons in Al crystal for protons (Zi = 1), as a function of the projectile velocity. Solid 
and open circles represent the results obtained in the RPA and the TDLDA, respectively. RPA 
and TDLDA stopping powers of electrons in an homogeneous electron gas with = 2.07 are 
represented by solid and dashed lines, respectively. 

FIG. 2. Contributions to the random stopping power of valence electrons in Al crystal coming 
from the excitation of (a) electron-hole pairs and (b) plasmons, versus the projectile velocity, as 
obtained from Eq. (3.14) within the RPA. Full band structure calculations are represented by solid 
circles. The corresponding jellium results are represented by solid lines. 

FIG. 3. Position-dependent stopping power of Al crystal for protons moving with v = 0.2, as 
obtained within the RPA. The projectile is incident in the (1 n m) direction and the stopping power 
is plotted as a function of the magnitude of the impact vector along the (a/? 7) direction, with 
(a) (lnm)=(100) and {a (3j)={0W); (b) (lnm)=(lll) and {a l3j)={-ll0). Solid lines represent 
full band structure calculations, as obtained from Eq. (3.16). Dashed curves represent LDA calcu- 
lations, obtained as explained in the text, and the horizontal dashed line represents the stopping 
power of a homogeneous electron gas with = 2.07. All results have been normalized to the value 
of the random stopping power of Al crystal for the same velocity. 

FIG. 4. Contour valence-density plots of the square lattice containing the projectile trajectory 
(Inm) and the impact vector {a(3j), with: (a) (lnm)=(100) and {a (3 j)={010); (b) (lnm)=(lll) 
and {a P j)={-l 10). Projectile trajectories with b = 1/4 are also represented. 

FIG. 5. Averaged valence-densities along the projectile paths depicted in Fig. 4, as a function 
of the impact parameter b. The solid line corresponds to a projectile trajectory along the (100) 
direction, as in Figs. 3a and 4a. The dashed line corresponds to a projectile trajectory along the 
(111) direction, as in Figs. 3b and 4b. All results have been normalized to the value of the average 
valence-density of Al, no- 
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Figure 2b 
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